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HORIZONTAL  AIR-MASS  TABLES 


ABSTRACT.  The  illumination  of  artificial  earth  satellites  dur¬ 
ing  their  motion  into  the  earth’s  shadow  is  a  problem  of  cur¬ 
rent  interest.  For  heights  above  the  surface  greater  than  20 
km,  the  tangential  solar  ray  incident  on  a  satellite  in  a  non- 
absorbed  region  of  the  spectrum  is  attenuated  principally  by 
Rayleigh  scattering.  The  attenuation  is  a  function  of  the 
total  number  of  molecules,  or,  if  the  molecular  weight  remains 
constant,  of  the  air  mass  along  the  path  of  the  ray  that  is 
tangential  to  any  level,  hQ,  above  the  surface. 

Tables  are  given  in  terms  of  mass  of  air  and  number  of 
molecules  for  a  vertical  path  originating  at  h  ,  a  horizontal 
path  tangent  to  h^,  the  horizontal  quantities  relative  to  a 
unit  atmosphere,  and  the  rate  of  change  with  respect  to  . 

The  tables  cover  an  altitude  range  of  0  to  100  km,  calculated 
for  an  atmosphere  terminated  at  200  km.  These  tables  are  in¬ 
tended  primarily  for  use  above  a  nominal  altitude  of  20  km, 
where  refraction,  which  was  not  allowed  for  in  the  computations, 
.  has  a  negligible  effect. 
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INTRODUCTION 

•  • 

• 

The  orbiting  of  large  balloon  reflectors  such  as  Echo  I  has  provided 
a  new  tool  for  investigating  the  composition  of  the  upper  atmosphere.  In 
addition  to  its  primary  utility  as  a  radio  communications  reflector,  the 
satellite  is  also  particularly  useful  as  an  extra-atmospheric^ reflector 
of  light  having  known  reflecting  properties.  Although  Echo  I  was  launched 
in  an  orbit  which  was  continuously  illuminated  by  the  sun,  the  orbit  has 
a  regression  rate,  relative  to  the  shadow  axis  of  the  earth,  of  approxi- 
,  mately  four  degrees  per  day.  Examination  of  the  geometry  involved  shows 
that  the^  orbital  plane  intersects  the  shadow  cone  of  the  earth  about 
60  percent  of  the  time .  On  these  ' shadow  passes '  the  solar  ray  illumin¬ 
ating  the  satellite  scans  the  atmosphere  just  before  the  satellite  enters 
the  shadow  and  just  after  it  leaves  the  shadow  on  each  rotation.  At  any 
time  that  the  illuminating  ray  is  traversing  the  atmosphere,  the  illumi¬ 
nance  at  any  wavelength  is  a  function  of  the  distribution  of  any  con¬ 
stituents  that  absorb  at  that  particular  wavelength,  and  is  also  a  func¬ 
tion  of  the  molecular  density  distribution,  which  alters  the  intensity 
by  Rayleigh  scattering.  Since  a  horizontal  path  through  any  atmospheric 
layer  is  much  longer  t'han  a  vertical  path,  it  is  possible  to  measure  the 
distribution  of  minor  atmospheric  constituents.  . 

Two  methods  of  analysis  are  available.  An  integral  technique  (Ref.  l) 
has  been  used  for  lunar  eclipse  measurements  of  the  ozone  distribution, 
and  has  been  modified  for  application  to  satellite  measurements  (Refs.  2 
and  3) •  Secondly,  a  matrix  method,  which  can  more  readily  be  programmed 
for  computers, 'has  also  been  devised  (Ref.  4). 

Although  the  methods  can  be  applied  to  the  determination  of  the 
atmospheric  density  distribution,  it  is  of  immediate  interest  to  calcu¬ 
late  the  Rayleigh  attenuation  in  ordpr  to  separate  the  absorption  and 
scattering  effects.  If  no  absorption  occurs  the  illuminance  at  the 
satellite  is  given  by 


i'U)  *•  i0(*)e 

where  Iq(x)  in  the  extraterrestrial,  solar  intensity  at  wavelength  X, 

k  (x)  is  the  Rayleigh  scattering  cross  section  or  coefficient,  and  M  is 
the  total  number  of  molecules  or  the  uir  mass  along  the  ray  from  the  sun 
to  the  satellite.  An  older  tabulation  of  M-vs-altitude  (Ref.  !>)  proved 
unsatisfactory  because  only  the  lower  levels  of  the  atmosphere  were  con¬ 
sidered.  Another  tabulation  (Ref.  0)  gives  air  muss  data  to  1^0  km  with 
rocket  density  data  obtained  at  White  Sands  up  to  19132.  Since  more  recent 
density  data  are  available,  we  felt  that  a  new  cal  culation  of  M  would  be 
useful . 
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It  vas  decided  to  neglect  the  effect  of.  refraction  in  the  tables 
A  included  in  this  report  because  in  the  region  of  interest  (20  to  100  km) 
the  errors  are  small.  However,  for  completeness  the  data  have  been  pre¬ 
sented  from  the  surface  to  100  km.  The  user  should  thus  exercise  some 
caution  in  use  of  the  data  at  the  lowest  levels.  We  hope  to  present  the 
data  corrected  for  refraction  at  a  later  date. 


CALCULATIONS 

We  can- define  'air  mass'  as  the  total  mass -per -unit  area  along  a 
specified  path  through  the  earth's  atmosphere.  The  particular  air  mass 
calculated  here  we  have  called  'horizontal  air  mass'.  The  path  in  this 
case  is  a  line  through  the  atmosphere  which  passes  the  earth  (assumed 
spherical)  at  some  minimum  distance,  h  ,  as  shown  in  Fig.  1.  . 


FIG.  1.  Horizontal  Air-Mass  Geometry. 
The  air  mass  M(h  )  is  then  given  by 


M(h  )  =  P  p(h)  dx 

u  '-Loo 


(1) 


where  p(h)  is  the  mass  density  of  the  air  at  height  h.  We  have  assumed 
that  thp  atmosphere  is  radially  symmetric.  Equation  (l)  can  then  be 
written  as 


because  of  this  symmetry. 


M(h  )  =  2 P  „p(h)  dx 
J  0 


(2) 
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T 


Of  X 


W«  now  wish  to  Introduce  h  as  the  variable  of  integration  in  place 
.  From  the  geometry  of  Fig.  1  ve  have  , 


x2  +  (R-H»0)2  -  (R+h)2 


(3) 


and  then 


(R+h)dh  _  (R4h)dh _  m 

-  (E4h af  V(h-h0)(ffl^0) 


Equation  (2)  now  reads 

M(hQ)  -2  f  — -  (5) 

ho  [ ( h-hQ) ( 2R+h4hQ) 

Although  Eq.  (5)  is  a  convergent  Integral  for  any  reasonable  func¬ 
tion  p(h),  it  presents  difficulties  if  the  integration  is  performed 
numerically  as  is  done  here.  The  upper  limit  presents  no  problem  since 
the  earth's  atmosphere  can  be  considered  to  terminate,  for  practical 
calculations,  at  some  finite  h. 

We  can  dispose  of  the  discontinuity  at  h  ■  hQ  by  the  following  sub¬ 
stitution: 


2 

y  - 


h-ho 


(6) 


Then 


dh  «  2ydy 


and  Eq.  (5)  becomes 


M(ho) 


p(y2+hQ) [R+hQ+y2]dy 
[2(R+hQ)+y2]^ 


(7) 
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For  practical  calculations  the  range  of  y  can  he  taken  as  # 


w  * y  * 0 


(8) 


where 


y  ■  [h  «h  ]4 
Jmax  max  o 


(9) 


We  shall  say  more  about  h  later. 

•  max 

Although  the  density  of  the  atmosphere  is  approximately  an  exponen-  . 
tial  function  of  height,  this  was  not  considered  a  sufficiently  accurate 
approximation  for  the  present  purpose.  Instead,  the  ARDC  1959  model 
atmosphere  (Ref.  6)  was  used.  The  properties  of  the  atmosphere  are 
tabulated  functions  of  height.  With  this  density  function  it  is  nec¬ 
essary  to  integrate  Eq.  (7)  by  numerical  methods. 

The  integration  was  performed  with  Slnqpson's  rule.  The  density 
data  of  the  ARDC  atmosphere  is  tabulated  at  various  intervals,  which  are 
larger  at  higher  altitudes.  This  fact,  and  the  fact  that  Simpson's  rule 
requires  equal  intervals,  necessitated  interpolation  for  the  values  of 
2 

p(y  +hQ).  The  equal  intervals  must  of  course  be  in  the  variable  of 

integration  y.  Second  order  interpolation  was  used.  The  numerical  work 
was  done  on  the  IBM  709 • 

The  upper  limit  for  h  was  set  at  either  200 

became  10  ^  of  the  initial  density  corresponding 

red  first.  The  atmospheric  density  was  taken  to 
The  range  of  hQ  was 


km  or  where  the  density 

to  h  ,  whichever  occur- 
o 

be  zero  above  200  km. 


0  S  hQ  £  100  km  (10) 


and  the  air  mass  was  computed  for  every  1,000  m  in  this  interval. 

Another  quantity  of  interest  in  these  calculations  is  the  'zenith 

air  mass'.  This  is  the  air  mass  computed  along  the  vertical  from  h.  to 

o 

200  km.  Thus 


4 


• 

• 
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• 

1 

• 

200,000 

p(h)dh 


(11) 


Again  the  integration  was  performed  by  Simpson's  rule.  In  this  integra¬ 
tion  no  interpolation  was  used.  The  interval  of  integration  was  the 

interval  of  tabulation  of  h  .  The  subroutine  which  performed  the  inte- 

o 

gration  took  care  of  the  end-point  automatically  when  there  was  an  even 
number  of  density  points  at  some  interval  of  hQ.  The  basic  Simpson's 

rule  requires  an  odd  number  of  points. 

The  'relative  air  mass'  was  also  computed.  We  define  the  relative 
air  mass  by 

M(h  ) 

H(ho)  ■  jg§>  (1S> 


that  is.  the  zenith  air  mass  at  h  ■  0  is  used  as  the  reference. 

•  o 

The  final  quantity  which  was  computed  was  the  derivative  of  the 
relative  air  mass  with  respect  to  height  hQ.  This  was  oomputed  using 

five-point  formulas  (Ref.  7). 

The  numerical  work  was  checked  by  several  means.  First,  the  zenith 
air  mass  for  hQ  should  be  equal  to  the  atmospheric  pressure  at  that  hQ 

with  suitable  unit  conversion,.  A  small  table  of  air  pressure  and  accel¬ 
eration  of  gravity  from  the  ARDC  atmosphere,  and  the  zenith  air  mass, 
calculated  pressure,  and  percentage  difference  between  the  pressures  are 
presented  in  Table  1.  The  error  of  about  2  in  1,000  does  not  seem  un¬ 
reasonable,  and  may  be  ascribed  to  errors  in  Simpson's  rule  integration 
and  other  numerical  processes.  The  large  error  at  hQ  «=  100  km  is  proba¬ 
bly  due  to  the  termination  of  the  atmosphere  at  200  km. 

A3  a  check  on  the  integration  for  the  horizontal  air  mass,  an  ex¬ 
ponential  density  was  assumed  which  passed  through  the  ARDC  density  at 
h  =  0  and  h  =  100  km.  If  we  set 

p(h)  =  Ae_ah,  (13) 
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* 


fEABLE  1.  Zenith  Air-Mass  Check 


(meters) 

him 

w 

(gm/cm2) 

calc 

(mh) 

Pressure 

difference 

0 

1.01325  x  103 

9.8067 

I.O356  x  103 

I.OI55  x  103 

+0.23$* 

20,000 

5.5293  x  101 

9.7452 

5.6851  x  101 

5.540  x  10+1 

+0.20$ 

50,000 

8.7858  x  10 

9.6542 

9.1223  x  10 ”1 

8.807  x  10 -1* 

+0.24$ 

70,000 

6.0209  x  10  "2 

9.59^2 

6.2864  x  10“2 

6.031  x  10"2 

+0.17$ 

90,000 

1.353  x  10"3 

9.535 

1.4200  x  10“3 

1.354  x  10"3 

+0.07$ 

100,000 

-4 

2.138  x  10 

9.505 

2.238  x  10 "4 

-4 

2.127  x  10 

-0.51$ 

Equation  (5)  becomes 


M(hQ)  =  2  /  — -^2 - 

ho  [ (h-hQ ) ( 2R+h+hQ ) F 


(14) 


n  J . 

where  A  -  1.2250  kg/nr  and  a  =  I.50  x  10 ”  /m.  Now  if  we  neglect  h  with 
respect  to  R, 


00 

M(ho)  syiit  Aj  2 - ^  (15) 

ho  (h-hQ)2 

In  n  strict  sense  this  approximation  may  not  seem  reasonable  since  we  are 
integrating  h  from  hQ  to  However,  'infinity'  here  is  on  the  order *of 

200  km,  as  compared  to  the  radius  of  the  earth  R  at  some  6,000  km.  The 

exponential  function  decays  faster  than  the  h2  term  increases.  In  any 
event,  we  are  after  an  estimate  to  check* the  order  of  magnitude  of  the 
numerical  work.  Set 


2 


x 


h-h 

o 


(16) 
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Thert  Eq#  (15)  jreada 

•  •  , 

^  l  -ah 

M(h  )  -  2^/5*  Ae“aho/'  e”5*  fix  — !1  *(17) 

0  ^0  V* 

Using  the  same  density  function,  the  zenith  air  mass  is 

.  Mz(h0)  -  A  J  e_ahdh  «  |  e“eilo  (18) 

ho  .  . 

and  the*  unit  air  mass,  M„(o),  is 

* 

Mz(0)«|  (19) 

The  relative  air  mass  is  therefore 

•  • 

M(h  )  _ _  . 

e-an°  (20) 

A  few  numerical  values  are  given  in  Table  2. 


TABLE  2.  Estimates  of  Relative  Air  Mass 


h 

(meters) 

Relative  .air  mass 
Exponential  density 

Calculated 

0 

77.5 

70.3 

20,000 

3.86 

4.34 

50,000 

4.28  x  10'2 

6.13  x  10"2 

70,000 

2.13  x  lO-3 

4.95  x  10"3 

100,000 

2.37  x  10"5 

1.74  x  10"5 
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Exact  agreement  is  not  expected,  and  this  calculation  was  performed  to 
check  the  order  of  magnitude. 

* 

A  third  check  was  made,  using  an  atmosphere  of  constant  density. 

The  integration  with  a  finite  thickness  can  easily  be  performed  since 
the  density  is  constant.  Here 

M(ho)  *  2  J'  pax  -  2  Xxp  (21) 

From  the  geometry  in  Fig.  2, 


FIG.  2.  Geometry  for  Check  Calcu¬ 
lation  With  Constant  Density. 


it  is  easily  seen  that 

x2x  =  (R+l^)2  -  (R+ho)2  =  (h1-ho)(2R+h1+ho)  (22) 


In  this  case 

M(ho)  =  2p[(h1-hQ)(2R+h1+ho)]'2' 


The  unit  zenith  air  mass  M  (0)  Is 

z 


pdh  =  ph^ 


(23) 


(2k) 
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M(h  )  0  x 

«w  -  ET  [<VV(2R+V*o>12  (25) 

Table  3  shows  values  of  the  relative  air  mass  computed  by  the  machine  pro 
gram  and  by  hand  from  Eq.  (25). 


TABLE  3*  Relative  Air  Mass  for  Constant  Density 


h 

(meters) 

Relative 

air  mass 

.^-difference 

Machine 

Eq-  25 . 

0 

16.077 

16.084 

0.044 

20,000 

15.261 

15.271 

0.066 

50,000 

13.945 

13.956 

0.079 

70,000 

12.982 

13.003 

0.162 

100,000 

11.409 

11.417 

.  0.145 

The  last  check  was  made  on  the  derivative  of  the  relative  horizontal 
air  mass.  This  was  done  with  the  same  formulas  that  were  used  in  the  * 
machine .  Agreement  -was  noted  to  the  last  figure  in  the  machine  tabula¬ 
tion  of  the  derivative  for  several  values  of  h.  . 

Strictly  speaking,  the  air  mass,  M,  con  be  used  for  calculating 
Rayleigh  scattering  only  if  the  molecular  weight  is  constant .  If  the 
latter  is  not  constant  the  number  density  is  the  proper  density  to  use. 
Calculations  of  the  total  number  of  moleeules  per  square  centimeter 
(integrated  number  density)  for  horizontal  and  Zenith  paths  were  also 
made,  using  Equations  (7),  (ll),  and  (12),  and  substituting  number 
density  for  mass  density. 

The  horizontal  air  mass  and  associated  quantities  are  presented  in 
Table  4.  The  two-digit  number  following  the.  entries  indicates  the  power 
of  ten  by  which  the  number  ■should  be  multiplied.  Thus,  .12250-02  = 
.12250x10"^.  The  air-raass  entries  may  be  considered  accurate  to  about 
0.2%  on  the  average  if  a  single  figure  may  be  quotfecl  for  this  quantity. 
The  derivative  error  is  estimated  to  be  less  than  2%. 

The  integrated  number  density  and  associated  quantities  are  given 
in  Table  5,  The  relative  horizontal  air  mass  from  Table  4  and  the 
relative  integrated  number  density  from  Table  5  depart  by  (at  most) 

0.3%  at  the  highest  altitudes.  For  most  practical  work  either  table  may 
be  used. 
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